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ABSTRACT 

We derive a theoretical relation between Rblr, the size of the broad- 
emission-line region of active galactic nuclei, and the observed soft X-ray 
luminosity and spectrum. We show that in addition to the well known 
Rblr ~ L^^'^ scaling, Rblr should depend also on the soft X-ray spec- 
tral slope and derive the expected relation between Rblr and the X-ray 
luminosity and spectral index. Applying this relation to calculate a 
predicted BLR radius for ten AGN with reverberation data, we show 
that including the dependence on the spectrum improves the agreement 
between the calculated BLR radius and the radius independently deter- 
mined from reverberation mapping. Similarly we evaluate an expression 
for the line width, and show that including the dependence on the spec- 
trum significantly improves the agreement between the calculated BLR 
velocity dispersion and the observed FWHM of the H/3 line. The theo- 
retical expression for the line width also provids a physical explanation 
to the anti-correlation between the soft X-ray slope and the emission-line 
width observed in narrow-line Seyfert galaxies. 



Subject headings: galaxies: active — galaxies: nuclei — galaxies: 
Seyfert — quasars: emission-lines — X-rays: galaxies 
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1. INTRODUCTION 

Recent results from reverberation-mapping 
of the broad emission-line regions (BLR) in 
AGN indicate that the BLR distance from 
the central radiation source roughly scales as 
r cx L^/"^ (Koratkar & Gaskell 1991;Kaspi 
et.al. 1996; Peterson 1995). Our model pro- 
vides a physical explanation to this scaling, 
and predicts an additional testable depen- 
dence on the spectral shape. 

Most workers today agree that the line 
width in AGN is induced by Kcplcrian bulk 
motions. If this is the case, the BLR dis- 
tance is directly related to the line width, and 
studying such relations could provide more 
observational input. Such an input seems 
to be found in a new class of Narrow-Line 
X-ray Galaxies (NLXG), which are similar 
to Seyfert 1 galaxies except for having rela- 
tively narrow permitted lines (FWHM < 
2000km/scc) , and often stronger Fell emis- 
sion (e.g. Walter and Fink 1993). Many 
NLXG show steep soft X-ray spectra, ax > 

2, while ordinary (broad hne) AGN such as 
Seyfert I's and quasars usually have ~ 
0.5 — 1.5. Mixed samples (NLXG and or- 
dinary AGN) show a strong correlation be- 
tween the H/3 line width and the soft X-ray 
spectral index (Boiler, Brandt and Fink 1996; 
Wang, Brinkmann and Bergeron 1996, here- 
after WBB). The latter authors found the soft 
X-ray spectral index is correlated also with 
the FeII/H/3 ratio, and with the H/5 equiva- 
lent width. 

We derive simple analytic relations be- 
tween soft X-ray continuum spectrum (lumi- 
nosity and spectral index) and the BLR size 
and the width of the broad emission lines, 
which explain the observed correlation be- 
tween the X-ray slope and the line width (sec- 



tion 2). Applying the model to samples of 

AGN it gives a good agreement with the BLR 
sizes determined by reverberation mapping, 
and with the observed FWHM(H/?) - cf. Wan- 
del and Boiler 1997b ( section 3). Finally, in 
section 4 we suggest explanations of two ad- 
ditional empirical correlations (related to the 
H/3 equivalent width and the Fell emission) 
and discuss the validity of our basic assump- 
tions. 

2. THE MODEL 

While there is probably a range of physical 
conditions (ionization parameter and density 
etc.) in the broad Hbeta line-emitting gas 
in AGNs, the emission of each line is domi- 
nated by emission in a relatively narrow op- 
timum range of conditions (Baldwin et al. 
1995). We will therefore make the assump- 
tion that the conditions relevant for H/3 can 
be approximated by a single value. We show 
that within the observational constraints on 
the far UV and soft X-ray bands of AGN spec- 
tra, the BLR size and the line width predicted 
by the model are not very sensitive to the 
detailed shape of the spectrum in the EUV 
band. We then estimate the BLR size, and 
assuming Keplerian velocity dispersion we re- 
late it to the line width. The new element in 
our scheme is the use of the soft X-ray spec- 
tral index and luminosity to estimate the ion- 
izing spectrum. The spectral index enters in 
the model as follows: a softer (that is steeper) 
spectrum has a stronger ionizing power, and 
hence the BLR is formed at a larger distance 
from the central source, has a smaller veloc- 
ity dispersion and produces narrower emission 
lines. We parametrize the form of the ionizing 
EUV continuum in terms of the the break be- 
tween the UV and X-ray bands and the mea- 
sured soft X-ray slope a^- 
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2.1. The BLR radius 

We assume that the spatial extent of the 
BLR may be represented by a characteristic 
size - e.g. the radius at which the emission 
peaks. 

The physical conditions in the line-emitting 
gas are largely determined by the ionization 
parameter (the ratio of ionizing photon den- 
sity to the electron density, Ue, e.g. Netzer, 
1990) U = Q,„„/47ri?2cnewhcre the ioniz- 
ing photon flux (number of ionizing photons 
per unit time) is Q,on = /^i^(^)f , F{E) 
is the luminosity per unit energy, and Eq —1 
Rydberg=13.6 eV. Defining the ionizing lumi- 
nosity. Lion = Se'o F{E)dE and the mean en- 
ergy of an ionizing photon, Eion = Lion/Qion, 
the BLR radius (expressed in light-days) may 
be written as 
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(1) 

where nio = Ue/lO^^ cm ^, L^u = L^/IO'^'^ 
ergs/sec is the observed X-ray luminosity, 
e = Eion/Eo is the mean photon energy in 
Rydbergs, and = eLx/Lion = Lx/EoQion- 
Typical values in the gas emitting the high 
excitation broad lines are U ~ 0.1 — 1 and 
Ue ~ 10^° — lO^^cm"^ (e.g. Rees, Netzer 
and Ferland, 1989), so that Uriio ~ 0.1 - 10. 
For Lx we consider two practical cases: broad 
band luminosity, (for the ROSAT band) Lx — 
lo ikev F{E)dE and a monochromatic lumi- 
nosity defined by Lx = EF{E) at E=0.3keV. 
The two cases give quite different Rblr vs. 
ax curves, as shown in fig.l. 

2.2. The ionizing spectrum 

The ionizing flux is dominated by the EUV 
continuum in the 1-10 Rydberg regime, where 
most of the ionizing photons are emitted. 



Since the the continuum in this range can- 
not be observed directly, we try to estimate 
it by extrapolation from the nearest observ- 
able energy bands. The far UV spectrum 
has been observed beyond the Lyman limit, 
for about 100 quasars (Zheng et al, 1996), to 
wavelengths of 600A . and for a handful lu- 
minous, high rcdshift quasars to wavelengths 
of 350A. The average spectrum has a slope of 
q; ~ 1 in the 1000-2000A band, and below 
900A it steepens to a ~ 2. In the soft X-ray 
band the observed spectrum is often steeper 
than this (especially in NLXGs, see Hasinger 
et al. 1993; Walter and Fink, 1993, but com- 
pare Laor et al. 1997). This indicated there 
may be a break or a turnover at some inter- 
mediate energy Ei, between 10 and lOOeV (cf. 
Mathews and Ferland 1987). We assume that 
the soft X-ray spectrum can be extrapolated 
to lower energies down to some break energy 
El,, and below the break we take E'"^. For the 
hard X-rays [E > 2keV) we use the "univer- 
sal" E~'^''^ power law, observed in most AGN 
with a high energy cutoff at lOOkeV (since the 
ionizing photon flux is dominated by the con- 
tinuum near 1 Ryd, the model is insensitive to 
the assumed hard X-ray spectrum). In sum- 
mary, we approximate the ionizing continuum 
spectrum by a broken power law of the form 



F{E) DC < 



Eq < E < Ef,; 
E,<E< E^; (2) 



Eh < E < Ejnax'i 

where the break energy Eb > Eq is a. free 
parameter, E^ is the break energy at the 
hard X-ray band, and E^ax is the high en- 
ergy cutoff. In the calculations below we take 
Et = O.lkeV and Eh = 2keV. 

Fig. 1 shows how the calculated BLR size 
depends on the X-ray spectral index, on the 
shape of the EUV spectrum (parametrized by 
Eh) and on the normalization of the X-ray 
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luminosity (broad-band or monochromatic). 
While the dependence on is strong {Rblr 
increases by a factor of 100 when increases 
from to 4), the dependence on the EUV 
shape is quite weak; changing Ef, in the rele- 
vant range (13.6-lOOeV) changes Rblr by less 
than 30% (for «^ < 3). As expected, we find 
that Rblr is almost independent of the hard 
X-ray parameters - and E^ax- 

2.3. Line width 

Assuming that the line width is induced by 
Keplerian motion in the gravitational poten- 
tial of the central mass, the velocity disper- 
sion corresponding to the full width at half 
maximum of the emission lines is given by 
V ^ ^JGM/R where M is the mass of the cen- 
tral black hole and r the distance of the broad 
line region from the central source. Eq. (1) 
gives 

lonn fUn^V'' ( ^_\''' ^ i 
V fa 1900 — z km/s. 

(3) 

In order to relate the unknown mass to the 
observed luminosity we assume that the cen- 
tral mass approximately scales with the lu- 
minosity (Dibai 1981; Joly et al. 1985; Wan- 
del and Yahil 1985; Wandel and Mushotzky 
1986). In terms of the Eddington ratio these 
authors find for large AGN samples L/ Lem ~ 
0.01 — 0.1. Within this distribution, bright 
objects tend to have slightly larger L/M ra- 
tios than faint ones (cf. Koratkar and Gaskell 
1991) roughly L/M oc L^'^ or M oc L^'*. 
Combining this with the correlation between 
the optical and the X-ray luminosities oc 
^a75±o.o5 ^j^j..gg ;^ggg. Mushotzky and Wandel 

1989) gives 

M 7 X IO^L.44 {^^Y§r) ' 



and with eqs. (2) and (3) the hne-width 

may be related directly to the observed X-ray 
luminosity and spectral index: 

v{FWHM)p^bmQr]e]l^{a:,)L]l^ km/s (5) 

where ri = {Uniofl^{L/ LEdd/^-^^)~^'^ com- 
bines all the unknown parameters; for the 
Unw ~ 0.1 - 10 and L/Lem ~ 0.01 - 0.1 
stated above we have 77 ~ 0.2 — 2. 

3. COMPARISON WITH THE DATA 
3.1. BLR radius 

We have Compared the BLR size calcu- 
lated from the model with the distance from 
the central source obtained by reverberation 
mapping of the H/3 line, for a sample of 10 
AGN for which reverberation and X-ray data 
were available (fig 2 and table 1). As we dis- 
cuss below, the agreement actually improves 
by taking the X-ray spectral slope into ac- 
count. 

In the model calculations we have used 
the spectral index from the power-law fit to 
the ROSAT 0.1-2.4 keV band and E^^lRyd, 
(with the exception of NGC 3783, for which 
we have used = 1.5, see Walter and Fink, 
1992; note also that NGC 4151 has an ex- 
tended emission, Warwick et.al. 1995). The 
horizontal error bars represent the combined 
error in the luminosity and in the spectral in- 
dex. Allowing for the uncertainty in Unio ~ 
0.1 — 1 (represented by the dashed lines in 
fig. 2), the agreement is quite good: all the 
points lie well within these boundaries. In 
order to test the significance of our model, we 
have calculated the BLR radii without taking 
into account the spectral dependence (that is, 
assuming all objects have the same soft X-ray 
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spectral slope, which we set to the sample av- 
erage, a ^=1.45). We find that the difference 
between BLR size (calculated using only the 
luminosity scaling) and the reverberation size 
is significantly larger when the spectral de- 
pendence is not taken into account for three 
objects (NGC4151, PG0844 and NGC 3783) 
shown as gray dots in fig 2. For the other 
objects the difference is small in both calcu- 
lations. 

In order to compare our model to the em- 

1 /2 

pirical R LJp^ relation, we have calculated 
the BLR radius using eq. (1) with L = Lopt 
and e = L We denote these three models by 
R{L^,Oi^), R{L^), and R{Lopt). 

The statistical significance of each model 
is tested by the statistic, 

, ^ ^ [log(i?,(obs)) - log(i?,(mod))]^ 
^ Y a2(obs) + f7Kmod) 

The variance of each object is taken as the 
sum of the squared standard deviations of the 
observed radius (the error in the reverbera- 
tion determination) and the caknilated radius 
(the errors in the involved observables - L^, «^ 
or Lopt- The values obtained for the three 
models are x^{Lx,ax) = 0.88, x^iLx) = 
1.77, and x^{Lopt) — 1-56. The correspond- 
ing confidence levels p are given hy 1 — p — 
r(z//2, X^/2) where v is the number of degrees 
of freedom (here taken as 9). For the three 
models we get l-p=0.0043, 0.0091 and 0.0077 
respectively. 

3.2. Line width-spectral index correla- 
tion 

Equation (5) predicts an explicit relation 
between the velocity dispersion (associated 
with the line width), the X-ray luminosity 
and the spectral index, namely a surface in 
the — V — Lx space. For a fixed value 



of Lj. this gives a curve in the v — plane. 
Figure 3 shows such curves of FWHM vs.ccj, 
for several values of the luminosity, = 
lO'^^-lO*^^ ergs/sec (cf. Wandel and Boiler 
1997). Overplotted are the data points - 
FWHM(H/3) vs. for a sample of AGN (see 
below). The model seems to reproduce the 
distribution of the data very well, and in par- 
ticular it explains the observed anticorrelation 
between the H/3 line width and the soft X-ray 
spectral index. 

3.3. Predicted vs. observed line width 

In order to test the model prediction over 
the three dimensional a — v — space we 
compare the observed line width to the value 
calculated with eq. (5) using the measured X- 
ray spectral index and luminosity for a sample 
consisting of 33 ordinary AGN from Walter 
and Fink (1993) combined with 32 NLXGs 
(Boiler, Brandt and ,Fink 1996). As can 
be seen in fig. 4a, the agreement is very 
good, and most of the objects fall well within 
the uncertainty strip of log FWHM(obs)=log 
v(mod)±0.5. Fig. 4b shows the same com- 
parison when the spectral information is not 
taken into account (but rather its sample av- 
erage, ax — 1.74): the correlation is signifi- 
cantly weeker. The correlation coefficient be- 
tween the observed and calculated line widths 
is r = 0.533 and 0.316 respectively. The t-test 
statistic for zero regression 

t = r^{N - 2)/^{l - r2) gives for the two 
cases t — 3.60 and 2.40 respectively, for which 
the confidence levels for zero correlation are 
1 —p = 0.0004 and 0.01 respectively, and the 
dependent confidence level (which measures 
the significance of the dependence on spectral 
index, after taking out the dependence on the 
luminosity) is p = 1 - 0.0004/0.01 = 0.96. 
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3.4. Equivalent width 

One may try to understand also the other 
strong correlations found in the data in terms 
of our model. The H/3 equivalent width 
is anti-correlated with the FWHM (Gaskell 
1985, WBB) and with the soft X-ray spec- 
tral index; AGN with a steep soft X-ray spec- 
trum have a lower EW(H/5) than flat- spec- 
trum AGN. To see how our model can ex- 
plain this, we recall that the equivalent width 
measures the fraction of the continuum flux 
reprocessed and emitted in the hne. We have 
shown (fig. 1) that the BLR distance from 
the central source increases with the spectral 
index. If the emission-line clouds are not pres- 
sure confined, it is reasonable that their sizes 
do not change significantly with the BLR ra- 
dius, hence if clouds are conserved, the solid 
angle covered by the clouds decreases with ra- 
dius, which implies that the equivalent width 
decreases with increasing (steepening) spec- 
tral index, as observed. 

The author thanks Th. Boiler, A. Laor, 
D. Maoz and H. Netzer for Tiseful discussions, 
M.I. for help with the statistical analyses and 
M. Gaskell for the careful refereeing. 
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Figure Captions 

Fig. 1. — The calculated BLR radius vs. 
the soft X-ray spectral index a. Shown are 
two assumed EUV spectra: (a) Ei, = IRyd, 
so that aEuv = olx (solid curves), and (b) 
= lOOeV, aEuv = 2 (dashed curves). 
For each case we show two luminosity nor- 
malizations: monochromatic - Lx{0.3keV) — 
lO^^crgs/scc (the upper two curves, indi- 
cated by O.SkcV) and the ROSAT band nor- 
malization (see text) - I/-r(0.1 — 2AkeV) = 
lO^^ergs/sec (the lower two curves, indicated 
by 0.1-2.4keV). 

Fig. 2. — The BLR size determined from re- 
verberation mapping of the H/? line vs. the ra- 
dius calculated from the model (eq. (1)) using 
the observed X-ray luminosity and spectral 
index. The reverberation radii are from Kaspi 
et.al. (1996), and the X-ray data from Walter 
and Fink (1993) except NGC 4151 (Warwick 
et al. 1995) and PG 0026 and PG 0844 from 
WBB. The grey circles represent the radius 
calculated from the model without taking the 
spectral dependence into account (sec text). 
The solid line represents Rrev = Rmodei for 
U niQ — 1 while the dashed parallel hues above 
and below it correspond to Uuiq = 0.1 and 10 
respectively. 

This 2-column preprint was prepared with the AAS 
IMgX macros v4.0. 



Fig. 3. — Theoretical curves of the X-ray 
spectral index vs. v(BLR) for fixed values 
of Lx = lO^^ergs/scc (lower curve) to = 
10"^^ ergs/sec (upper curve) superimposed on 
the data for the sample of NLXGs and normal 
Seyferts. 

Fig. 4.— Observed FWHM(H/5) vs. 
v{Lx, Oix) (the velocity dispersion calculated 
with the observed X-ray spectral index and 
luminosity, eq. (5)) - figure 4a, and 
FWHM(H/5) vs. v{L.x) (calculated with only 
the observed X-ray luminosity- figure 4b. The 
diagonal line represents the equality {v — 
FWHM) for Tj = 0.6, while the dashed par- 
allel lines above and below the diagonal are 
for = 0.2 and 2 respectively. 
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Table 1 

X-RAY AND Optical data, reverberation and model BLR radius. 



Object 


Name 


Log Lx 


Log Lopt 












1 


NGC 3783 


43.5±0.16 


43.5±0.2 


1.50±0.40 


0.92±0.08 


0.91±0.29 


0.88±0.08 


0.90±0.1 


2 


NGC 5548 


44.0±0.06 


43.8±0.2 


1.21±0.15 


1.23±0.10 


l.OOiO.ll 


1.13±0.03 


1.05±0.1 


3 


MKN 279 


44.26±0.07 


44.2±0.1 


1.15±0.17 


l.OOibO.lO 


1.10±0.13 


1.26±0.035 


1.24±0.05 


4 


AK 120 


44.86±0.08 


44.15±0.25 


1.63±0.29 


1.46±0.15 


1.66±0.19 


1.56±0.04 


1.20±0.12 


5 


MKN 590 


44.2±0.08 


44.3±0.2 


1.01±0.24 


1.27±0.04 


1.00±0.17 


1.23±0.04 


1.30±0.1 


6 


PG 0804 


45.2±0.1 


45.3±0.1 


1.53±0.41 


1.92±0.23 


1.78±0.26 


1.73±0.05 


1.80±0.06 


7 


NGC 4541 


43.1±0.1 


43.2±0.2 


1.80±0.1 


0.85±0.25 


0.87±0.10 


0.68±0.05 


0.74±0.1 


8 


PG 0953 


45.16±0.05 


45.7±0.2 


1.31±0.28 


1.95±0.20 


1.64±0.17 


1.71±0.025 


1.98±0.1 


9 


PG 0844 


43.7±0.08 


45.4±0.3 


1.82±0.95 


1.57±0.27 


1.17±0.50 


0.98±0.04 


1.83±0.15 


10 


NGC 3516 


42.7±0.03 


43.5±0.1 


1.49±0.42 


0.59±0.12 


0.50±0.23 


0.48±0.015 


0.88±0.05 



"^Log of BLR radius from reverberation mapping in light days 

''Log of BLR radius calculated from the model in light days (see text) 



8 




log R (model) [light-days] 




2000 4000 6000 8000 

FWHM (km/s) 



-^ 2.8 
2 . 6 





no spectral 




dependence 


• * • • 
t 

' ' • * • 




t 

• • • 




• • • , 

••• • *: • 



2.6 2.8 3 3.2 3.4 3.6 3.8 4 



log V (model) 



